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Abstract We present a first detailed spectrum synthe- 
sis calculation of a supernova-fallback disk composed of 
iron. We assume a geometrically thin disk with a radial 
structure described by the classical a-disk model. The 
disk is represented by concentric rings radiating as plane- 
parallel slabs. The vertical structure and emission spec- 
trum of each ring is computed in a fully self-consistent 
manner by solving the structure equations simultane- 
ously with the radiation transfer equations under non- 
LTE conditions. We describe the properties of a specific 
disk model and discuss various effects on the emergent 
UV/optical spectrum. 

We find that strong iron-line blanketing causes broad 
absorption features over the whole spectral range. Limb 
darkening changes the spectral distribution up to a fac- 
tor of four depending on the inclination angle. Conse- 
quently, such differences also occur between a blackbody 
spectrum and our model. The overall spectral shape is 
independent of the exact chemical composition as long 
as iron is the dominant species. A pure iron composi- 
tion cannot be distinguished from silicon-burning ash. 
Non-LTE effects are small and restricted to few spectral 
features. 

Keywords Radiative transfer; scattering • Neutron 
stars • hifall, accretion, and accretion disks 

PACS 95.30.Jx • 97.60.Jd • 98.35.Mp 



1 Introduction 

Anomalous X-ray pulsars (AXPs) are slowly rotating 
(Pj.Qt = 5 — 12 s) young (< 100 000 yr) isolated neutron 
stars. Their X-ray luminosities (« 10^^ erg/s) greatly ex- 
ceed the rates of rotational energy loss (« 10"^^ erg/s). It 
is now generally believed that AXPs are magnetars with 
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magnetic field strengths greater than lO^"' G and that 
their X-ray luminosity is powered by magnetic energy 
(Woods & Thompson 2006). As an alternative explana- 
tion the X-ray emission was attributed to accretion from 
a disk that is made up of supernova-fallback material 
(van Paradijs etal. 1995, Chatterjee etal. 2000, Alpar 
2001). The fallback-disk model has difficulties to explain 
IR/optical emission properties of AXPs. When compared 
with disk models, the faint IR/optical flux suggests that 
any disk around AXPs must be very compact (e.g. Perna 
etal. 2000, Israel etal. 2004). 

The discovery of optical pulsations in 4U 0142+61 
which have the same period like the X-ray pulsations 
(Kern & Martin 2002) appears to be a strong argument 
against the disk model. It was argued that reprocessing of 
the pulsed NS X-ray emission in the disk cannot explain 
the high optical pulsed fraction, because disk radiation 
would be dominated by viscous dissipation and not by 
reprocessed NS irradiation (Kern & Martin 2002) . Ertan 
& Cheng (2004), on the other hand, showed that these 
optical pulsations can be explained either by the magne- 
tar outer gap model or by the disk-star dynamo model. 
Therefore, the observation of optical pulsations is not an 
argument against the disk model. A spectral break in 
the optical spectrum of 4U 0142-1-61 was discovered by 
HuUeman etal. (2004). This was also taken as an argu- 
ment against the disk model because the authors do not 
expect such strong features from a thermally emitting 
disk. The recent discovery of mid-IR emission from this 
AXP (Wang et al. 2006), however, has strongly rekindled 
the interest in studies of fallback-disk emission proper- 
ties. While this mid-IR emission was attributed to a cool, 
passive (X-ray irradiated) dust debris disk by Wang et al. 
(2006), it was shown by Ertan on this conference that it 
can be explained with a model for an active, dissipating 
gas disk. If true, then the disk emission properties allow 
to conclude on important quantities, e.g., the magnetic 
field strength of the neutron star can be derived from the 
inner disk radius. 

Independent hints for the possible existence of fall- 
back disks come from pulsars with particular spin-down 
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Table 1 Characteristics of the nine rings that compose the disk model. Surface mass density S and emergent flux (expressed 
as Tefi) follow from the a disk prescription. The Rosseland optical depth thoss at the disk midplane, the mass density p and 
gravity log g at optical depth unity follow from our computations of the detailed vertical ring structure. The last column 
denotes the fraction of the disk area that is made up by each ring model in order to compute disk-integrated spectra. 



Ring 


R 

[1000 km] 


E 
[g/cm^l 


[1000 K] 




log p(rRoss = 1) 
[g/cm^] 


logff(TRoss = 1) 

[cm/s^] 


% area 
fraction 


1 


2.0 


2.9 


305 


3.1 


-6.1 


7.7 


0.0025 


2 


2.5 


2.8 


258 


3.2 


-6.2 


7.5 


0.01 


3 


3.5 


2.7 


201 


3.3 


-6.4 


7.2 


0.034 


4 


6.0 


2.6 


135 


3.5 


-6.0 


6.7 


0.084 


5 


9.0 


2.6 


100 


3.7 


-5.9 


6.4 


0.19 


6 


14 


2.5 


72 


3.3 


-6.4 


6.0 


0.62 


7 


25 


2.3 


46 


3.5 


-6.9 


5.5 


1.70 


8 


40 


2.2 


33 


3.9 


-7.1 


5.1 


33 


9 


200 


1.9 


9.8 


3.8 


-7.6 


3.6 


64 



properties. For example, the discrepancy between the 
characteristic age and the supernova age of the pulsar 
B1757-24 was explained by the combined action of mag- 
netic dipole radiation and accretion torques (Marsden 
etal. 2001). Even more, the presence of jets from pulsars 
such like the Crab and Vela can possibly be explained 
by disk-wind outflow interacting with and collimating 
the pulsar wind (Blackman & Perna 2004). 

A fallback-disk model was proposed in order to ex- 
plain the X-ray enhancement following a giant flare of 
the Soft Gamma Repeater SGR1900-hl4 (Ertan & Al- 
par 2003). The X-ray light curve is interpreted in terms 
of the relaxation of a fallback disk that has been pushed 
back by the gamma-ray flare. This model can also explain 
the long-term X-ray and IR enhancement light curves of 
the AXP IE 2259-1-58 following a major bursting epoch 
(Ertan etal. 2006). 

The presence of a fallback disk around the stellar 
r(^mnant of SN 1987A has been invoked in order to ex- 
plain its observed lightcurvc which deviates from the the- 
oretical one for pure radioactive decay (Mcyer-Hofmeister 
1992). From the non-detection of any UV/optical point 
source in the supernova remnant, however, tight con- 
straints for the disk extension can be derived (Graves 
etal. 2005). 

To our best knowledge, the emission from fallback 
disks in all studies was hitherto modeled with blackbody 
spectra. In view of the importance of disk models for 
the quantitative interpretation of observational data it 
is highly desirable to construct more realistic models by 
detailed radiation-transfer calculations. 



2 Radial disk structure 

For the modeling we employ our computer code AcDc 
(Nagel etal. 2004), that calculates disk spectra under 
the following assumptions. The radial disk structure is 
calculated assuming a stationary, Keplerian, geometri- 
cally thin a-disk (Shakura & Sunyaev 1973). As pointed 
out by Menou etal. (2001), for a comparison with ob- 



servational data one probably has to use a more elabo- 
rate model, because near the outer disk edge the viscous 
dissipation and hence the surface mass density decline 
stronger with increasing radius than in an a-disk. How- 
ever, the purpose of the present paper is to look for differ- 
ential effects of various assumptions. Qualitatively, these 
effects can be expected to be independent of the detailed 
radial disk structure. In any case, it would be no prob- 
lem to carry out the computations presented here with 
different radial structures. 

The a-disk model is fixed by four global input pa- 
rameters: Stellar mass M^, and radius i?* of the accretor, 
mass accretion rate Af, and the viscosity parameter a. 
For numerical treatment the disk is divided into a num- 
ber of concentric rings. For each ring with radius R our 
code calculates the detailed vertical structure, assuming 
a plane-parallel radiating slab. 

In contrast to a (planar) stellar atmosphere, which 
is characterized by Tgs and log 5, a particular disk ring 
with radius R is characterized by the following two pa- 
rameters, which follow from the global disk parameters 
introduced above. The first parameter measures the dis- 
sipated and then radiated energy. It can be expressed in 
terms of an effective temperature Teg: 

T^siR) = [1 - 3GM^M/8aTTR^. 

The second parameter is the half surface mass density E 
of the disk ring: 

S{R) = [1 - (i?^/i?)l/2] M/3TTW. 

a and G arc the Stefan-Boltzmann and gravitational 
constants, respectively, w is the depth mean of viscosity 
w{z), where z is the height above the disk mid-plane. The 
viscosity is given by the standard a-parametrization as 
a function of the total (i.e. gas plus radiation) pressure, 
but numerous other modified versions are used in the 
literature. We use a formulation involving the Reynolds 
number Re, as proposed by Kriz & Hubeny (1986). We 
chose Re = 15 000 which corresponds to a « 0.01. 

For the results presented here we selected the fol- 
lowing input parameter values. The neutron star mass 
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2000 4000 6000 8000 

R/R. 

Fig. 1 Radial disk structure: Effective temperature (top 
panel) and Keplerian rotation velocity (bottom panel). Ar- 
rows mark a reference point at a distance of R—40 000 km 
(about 4000 stellar radii R*) from the neutron star as dis- 
cussed in the text. 



is 1.4 M0. The radii of the inner and outer disk edges 
are 2000 and 200 000 km, respectively. The disk is rep- 
resented by nine rings or, more precisely, by nine radial 
grid points. The radiation integrated over the whole disk 
is then computed by assigning a weight to each point's 
spectrum that resembles the area fraction that it repre- 
sents. The main characteristics of the disk at the radial 
grid points are given in Tab. ^ The mass-accretion rate 
was set to M = 3 • 10^^ Mq/ji. Fig. ^ shows the radial 
run of Teff. We also display the Keplerian rotation ve- 
locity for the later discussion of our results. The radial 
distance from the neutron star is expressed in units of 
the NS radius which is set to i?*=9.7 km. But note from 
the above equations that the disk model is essentially 
independent of the stellar radius for large distances from 
the neutron star. 

While S{R) and roft(-R) in columns 3 and 4 of Tab.[T] 
follow from the a-disk assumption, the quantities in the 
next three columns are the result from our detailed verti- 
cal structure calculations described below. It shows that 
the entire disk model is optically thick. The Rosseland 
optical depth at the disk midplane r™^^'''"" is >1000 
at all radii. We also tabulate the mass density and the 
gravity at unity optical depth. That demonstrates that 
the conditions in the line forming regions of the disk re- 
semble those in white dwarfs at the inner disk radii up to 
main sequence stars at the outer disk radii. The strength 
of Stark line broadening therefore strongly depends on 
the distance of the emitting region from the neutron star. 



3 Vertical disk structure 

The vertical structure of each disk ring is determined 
from the simultaneous solution of the radiation transfer 
equations plus the structure equations. The latter ones 
invoke radiative and hydrostatic equilibrium plus charge 
conservation. The structure equations also consist of the 
non-LTE rate equations for the atomic population densi- 
ties. The solution of this set of highly non-linear integro- 
differential equations is performed using the Accelerated 
Lambda Iteration (ALI) technique (Werner & Husfeld 
1985, Werner 1986, Werner etal. 2003). 

The total observed disk spectrum, which depends on 
the inclination angle, is finally obtained by intensity in- 
tegration over all rings accounting for rotational Doppler 
effects. 



3.1 Radiation transfer, hydrostatic and radiative 
equilibrium 

We consider the radiation transfer equation for the in- 
tensity 1 1, at frequency v. 

A*— ^^^^ = -Ki,{z)Ii,{fi,z) + ri^{z) 

with the opacity Ki, and the emissivity rj^. z measures 
the geometrical height above the disk midplane and [l 
is the cosine of the inclination angle i. The equation is 
solved using a short characteristics method. Opacities 
and emissivities are computed using atomic population 
densities that are obtained by solving the non-LTE rate 
equations. Our code allows for the irradiation of the disk 
by the central source, however, the results presented here 
are computed with zero incident intensity. 

The radiation-transfer equations plus vertical struc- 
ture equations are solved like in the stellar atmosphere 
case, but accounting for two basic differences. First, the 
gravity (entering the hydrostatic equation for the total, 
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Table 2 Summary of non-LTE model atoms for silicon, sul- 
fur, and iron. The numbers in brackets at the iron ions give 
the number of individual lines summed up into superlines. 
Employed for a specific test run, the silicon and sulfur model 
atoms are tailored to the conditions encountered in disk ring 
number 8. 



T [1000 K] 



element 


ion 


NLTE levels 


lines 




Si 


III 


6 


4 






IV 


16 


44 






V 


1 







S 


III 


1 









IV 


6 


4 






V 


14 


16 






VI 


1 







Fe 


I 


7 


25 


(141821) 




II 


7 


25 


(218 490) 




III 


7 


25 


(301981) 




IV 


7 


25 


(1027 793) 




V 


7 


25 


(793 718) 




VI 


8 


33 


(340 132) 




VII 


7 


22 


(86 504) 




VIII 


7 


27 


(8 724) 




IX 


7 


25 


(36 843) 




X 


7 


28 


(45 229) 




XI 


1 










5 


10.2 


20.8 


42.5 


86.7 


176.6 


350.6 


734.5 


1500 



E 

i 1500 - 




2.0x10= 



Fig. 3 This cut perpendicular to the midplane shows the 
temperature structure of the disk. Note that the vertical 
scale (height above the midplane) is expanded. The height- 
to-radius ratio is about 0.015. The cross marks the depth at 
R=40 000 km where 



3.2 Non-LTE rate equations 



i.e. gas plus radiation, pressure) is not constant with 
depth, but increases with z. The gravity is the vertical 
component of the gravitational acceleration exerted by 
the central object (self-gravitation of the disk is negligi- 
ble): 

Second, the energy equation for radiative equilibrium 
balances the dissipated mechanical energy and the net 
radiative losses: 

POO 

9/4 pw GE/R^ = 47r / {rj^ - K^J„)dv, 
Jo 

where p and J^, denote mass density and mean inten- 
sity, respectively. In the case of a stellar atmosphere the 
left-hand side of this equation vanishes and we get the 
usual radiative equilibrium equation. The solution is ob- 
tained by a generalized Unsold-Lucy scheme and yields 
the vertical temperature structure. 

Having calculated the vertical structures and spectra 
of the individual disk rings, the ring spectra are inte- 
grated to get the spectrum of the whole accretion disk: 

i?,o 27r 

= cos(i) J J Iv{i^4)^r)r d(j)dr . 
Ri 

Here, Ri and Ro denote the inner and outer radius of 
the disk, and is the azimuthal angle. At this stage, the 
Keplerian rotation velocity Vrot is taken into account by 
assigning a Doppler shift of Av = ^ymt sini/isinz to the 
intensity emerging from a specific azimuthal ring sector. 



For each atomic level i the rate equation describes the 
equilibrium of rates into and rates out of this level and, 
thus, determine the occupation numbers nf. 

The rate coefficients Pij have radiative and electron coUi- 
sional components: Pij — Rij + Cij . The radiative down- 
ward rate for example is given by: 

(Jijiv) is the photon cross-section and (rii/nj)* is the 
Boltzmann LTE population ratio. 

The blanketing by millions of lines from iron aris- 
ing from transitions between some 10^ levels can only 
be attacked with the help of statistical methods (Ander- 
son 1989, Dreizler & Werner 1993). At the outset, model 
atoms are constructed by combining many thousands of 
levels into a relatively small number of super levels. The 
respective line transitions are grouped into superlines 
connecting these superlevels. In this case, the popula- 
tion densities of the superlevels are computed from the 
rate equations, in which the photon cross-sections crij(v) 
in the radiative rates Rij do not contain only a single line 
profile but all individual lines that are combined into a 
superline. As an example we show such a cross-section in 
Fig.[21 The complete spectrum of our disk model (A=4- 
300000 A) is sampled by 30 700 frequency points. 

The model atoms that we have created for our disk 
calculations are summarized in Tab.|21 Most important 
is the iron model atom. It comprises the first eleven ion- 
isation stages and a total number of more than 3 million 
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1000 2000 3000 4000 5000 6000 
wavelength / A 



Fig. 4 Relative contribution of single disk rings to the to- 
tal disk flux. For clarity, we only show the fluxes from rings 
number 1, 6, 8, and 9). The flux from ring 8 at i?=40 000 km 
with refi=33 000 K dominates the total disk spectrum at 
UV/optical wavelengths. 



lines. Atomic data are taken from Kurucz (1991) and the 
Opacity and Iron Projects (TIPTOPbase^). 



3.3 Disk composition 

The chemical composition of the supernova fallback ma- 
terial in the disk is not exactly known. It depends on 
the amount of mass that goes into the disk. A disk with 
a small mass (say < 0.001 Mq) will be composed of 
silicon-burning ash (Menou etal. 2001). For simplicity, 
the results presented here are obtained by assuming a 
pure-iron composition. It turns out that the emergent 
spectrum is insensitive against the exact composition as 
long as iron is the dominant species fSect. lTT|l . For a re- 
spective test run for one specific ring we assumed a com- 
position that represents a silicon-burning ash. It contains 
iron (80% mass fraction) as well as silicon and sulfur by 
10% each. 



3.4 Disk model properties 

Fig. 13 displays the temperature structure of the disk. 
The temperature varies between 1.5 million K at the mid- 
plane at the inner disk edge down to 6000 K in the upper 
layers at the outer disk edge. At all radii the vertical run 
of the temperature decreases almost monotonously with 
height above the midplane. A mild temperature reversal 
in the uppermost disk layers occurs. This turns out to be 
a non-LTE effect, because the respective LTE model has 
a strictly monotonous temperature run. We will discuss 
the consequences of this effect in Sect. 14.21 

Which disk regions contribute to the total disk spec- 
trum and to what extent? In Fig. 0] we plot the emergent 

^ |http://vizier.u-strasbg.fr/topbase/| 
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Fig. 5 Vertical stratification of the disk at i?=40 000 km 
(ring 8, Tcg^SS 000 K). The quantities are plotted against the 
column mass density that is measured from the disk surface to 
the midplane (from left to right). The vertical line at log m = 
—0.9 denotes the depth where tross = 1. The dashed curve in 
the top panel shows the temperature structure of a respective 
LTE model. 



astrophysical flux from the area of four disk rings (rings 
1, 6, 8, and 9, see Tab. P), i.e., the computed flux per 
cm^ is weighted with the ring area. The spectral flux 
distribution of the innermost ring with Tofi=305 000 K 
has its peak value in the soft X-ray region. The contri- 
bution of this innermost region to the optical/UV spec- 
trum is negligible. Cutting of the disk at this inner ra- 
dius (i?=2000 km), therefore, is justified if this spectral 
range is of interest. The disk region that is dominat- 
ing the UV/optical flux is represented by ring 8 with 
Toff =33 000 K. Its radius is 40 000 km, that is about 
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Fig. 6 Vertical iron ionisation structure at _R=40 000 km 
(ring 8, Tofi^SSOOO K). The vertical line at log m = -0.9 
denotes the depth where tross = 1. The dominant iron opac- 
ities in the line forming regions are from Fe III- VI. 
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Fig. 7 Detail from the emergent disk spectrum at 
ii=40 000 km (ring 8, Toff =33 000 K). We compare a pure 
iron composition (black line) with a silicon-burning ash com- 
position (red line). The differences are marginal. 



4000 neutron star radii. Its spectrum is dominated by 
strong blends of the numerous iron lines. Further out in 
the disk the effective temperature decreases and the flux 
contribution to the UV/optical spectrum decreases, too. 
Our outermost ring (number 9) has Tcff=9800 K, its flux 
maximum is at A=4000 A and it is fainter than the inner 
neighbor ring 8 over the whole spectral range. Cutting off 
the disk at this outer radius (i?=200 000 km) therefore 
does not affect the UV/optical spectral region. 



4 Results 

Because ring 8 is a representative disk region that deter- 
mines the UV/optical spectrum, we will discuss its prop- 
erties in more detail in the following three subsections, 
addressing the effects of chemical composition, non-LTE 
physics, and limb darkening on the spectrum. The radius 
of this ring is marked by arrows in Fig. ^ and by a cross 
in Fig. El 

In Fig. |S1 we show the vertical structure of the disk 
at this radius (i?=40 000 km). We plot the run of sev- 
eral quantities on a column-mass scale, measured inward 
from the surface toward the midplane of the disk. The 
vertical line at logm = — 0.9 marks the depth at which 
tross = 1, i-e., the region of spectral line formation. The 
temperature shows a strong increase towards the mid- 
plane and, as already mentioned, a mild temperature re- 
versal in the optically thin surface layers. Together with 
the temperature run, the panels showing the mass den- 
sity and gravity structure indicate that their values in 
the line forming region are comparable to those encoun- 
tered in hot subdwarfs. The lowest panel shows that the 
disk height H at this distance from the NS is « 600 km, 
i.e. RjH = 0.015. 

In Fig. |H| we show the vertical ionisation structure 
of iron in the disk at i?=40 000 km. The dominant ion- 
isation stages in the line forming regions are Felll-VI. 



At the midplane FeVII is dominant. The temperature 
at any depth is so high that Fel-II do not signiflcantly 
contribute to the spectrum. 

4.1 Effects of chemical composition 

In Fig.[3we show the flux spectrum of ring 8 in the wave- 
length range A=6200-6500 A. It has been calculated for a 
pure iron composition as well as for a Fe/Si/S=80/10/10 
composition representing silicon-burning ash. The differ- 
ence between the spectra is very small, because they are 
completely dominated by the extremely large number of 
iron lines. We conclude that the exact disk composition 
is not affecting the spectrum as long as iron is the dom- 
inant species. 

While silicon and sulfur do not affect the overall spec- 
trum by continuous background opacities, line features 
can be seen in the computed spectra, e.g. the Si IV reso- 
nance line in the UV. The line depth reaches about 50% 
of the continuum level but it would be difficult to detect 
even in medium resolution spectra when the disk incli- 
nation is high and the spectra are Doppler broadened by 
rotation. 



4.2 Significance of non-LTE effects 

For our particular disk model we expect that non-LTE ef- 
fects are not very large. This is because of the relatively 
high gravities in the line forming regions, ranging be- 
tween log 3=3.6 in the outermost ring with rcff=9800 K 
and log 3=7.7 in the innermost region with 305 000 K. In 
Fig.|H|we compare the spectra of an LTE and a non-LTE 
model of disk ring 8 in the wavelength range A=1800- 
2900 A, where the largest deviations were found. Indeed, 
non-LTE physics affects only narrow spectral regions. 
Only there, fiux differences occur to an extent that the 
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2000 2200 2400 2600 2800 
wavelength / A 

Fig. 8 Comparison of LTE and NLTE emergent disk intensi- 
ties (red and blue lines, respectively) at i?=40 000 km (ring 8, 
i = 87°, Toff =33 000 K). 



equivalent width of line blends changes by a factor of two. 
Accordingly, the temperature structures of both models 
deviate only in the uppermost layers of the disk (see top 
panel of Fig. [SJ and, hence, only strong spectral lines 
that are still optically thick can be affected. 

4.3 Limb darkening effects 

Our model spectra show distinct limb-darkening effects. 
The situation is similar to the stellar atmosphere case 
(center-to-limb variation of the specific intensity). Look- 
ing face-on we see into deeper and hotter (and thus 
"brighter" ) layers of the disk when compared to a more 
edge-on view. In Fig. |^ we compare the specific inten- 
sity emitted by ring 8 (per unit area) for a high and a 
low inclination angle. Overall, the "edge-on" spectrum is 
roughly a factor of two fainter than the "face-on" spec- 
trum in the optical region. The difference increases to- 
wards the UV and amounts to a factor of about three. 
We conclude that limb darkening effects are important 
when disk dimensions are to be estimated from magni- 
tude measurements. 

It is also interesting to compare the intensities with 
a blackbody spectrum (Fig.lHl). Depending on the wave- 
length band, the blackbody over- or underestimates the 
"real" spectrum up to a factor of two in the optical and 
a factor of four in the UV. 



4.4 Rotational broadening 

We have seen that the spectrum of an iron-dominated 
disk is characterized by strong blends of a large number 
of lines. At some wavelengths broad spectral features ap- 
pear. It remains to be seen if Doppler effects from disk 
rotation smears out these features or if they could still 
be detectable. From the bottom panel of Fig. ^ we see 
that rotational broadening amounts to w sin i- 2000 km/s 




1000 2000 3000 4000 5000 6000 
wavelength / A 



Fig. 9 Effect of limb darkening: Specific intensity of the disk 
at 7?=40 000 km (ring 8) seen under two inclination angles, 
namely 87° (i.e. almost edge-on, blue) and 18° (i.e. almost 
face-on, red). For comparison we also show a blackbody spec- 
trum with T=Tcff=33 000 K. 

at i?=40 000 km, corresponding to an orbital period of 
about two minutes. When seen almost edge-on, this ro- 
tational velocity is equivalent to a Doppler broadening 
of about Z\A=25 A at A4000 A which clearly smears out 
any individual line profiles. In Fig. llOl we display the ro- 
tationally broadened spectrum of the entire disk model, 
seen under three different inclination angles. It is obvi- 
ous that the broad line blends are so prominent that they 
do not disappear even for an almost edge-on view of the 
disk. 

Among the strongest features is a 200 A wide line 
blend at A1500 A with an absorption depth of about 50% 
relative to the continuum. Should the disk be cooler, then 
disk regions with Tcff w90Q0 K could dominate the optical 
emission and the disk spectrum might look more like that 
emitted by ring 9 in our model (Fig.^. Strong iron-line 
blanketing could cause a spectral break. This contrasts 
with a statement in HuUeman etal. (2004), where the 
spectral break observed in the optical energy distribution 
of the AXP 4U 0142+61 is suggested as an argument 
against disk emission. 



5 Summary and outlook 

We have computed a model for a supernova-fallback disk 
in order to study its structure and optical/UV emission 
properties. We assumed an a-disk for the radial struc- 
ture and performed detailed non-LTE radiation transfer 
calculations for the vertical structure and spectrum syn- 
thesis. The input parameters were: 

Neutron star mass: 1.4 Mq 

Inner and outer disk edge radii: i?=2000 and 200 000 km 
Mass-accretion rate: M = 3 • 10^^ MQ/yr. 

We have identified that the disk region in the vicinity 
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Fig. 10 Complete disk spectrum including Kepler rotation 
broadening, seen under three inclination angles. The broad 
iron-line blends are detectable even in the almost edge-on 
case. 



star. This will reveal the relative importance of viscous 
dissipation and reprocessed irradiation that is discussed 
in the context of simultaneous optical and X-ray pulsa- 
tions in the AXP 4U 0142+61. At the moment we do not 
dare to make any prediction how this affects the results 
presented here. 

The innermost disk ring has a very high effective tem- 
perature and its flux distribution peaks in the soft X-ray 
band. It needs to be investigated systematically under 
which conditions (inclination, inner disk radius, accre- 
tion rate) the innermost disk regions can contribute to 
the thermal spectrum of the magnetars. 
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of i?=40 000 km is the main contributor to the total disk 
spectrum at UV/optical wavelengths. We therefore in- 
vestigated in some detail the disk properties at this ra- 
dius. 

We summarize our results as follows: 

— The overall disk spectrum is independent of the de- 
tailed chemical composition as long as iron is the 
dominant species. In particular, a pure-iron compo- 
sition is spectroscopically indistinguishable from a 
silicon-burning ash composition. 

— The overall disk spectrum is hardly influenced by 
non-LTE effects, however, equivalent widths of indi- 
vidual line blends can change by a factor of two. 

— Limb darkening affects the overall disk spectrum (in 
addition to the geometric foreshortening factor cos i). 
Depending on inclination and spectral band, the disk 
intensity varies up to a factor of three. 

— Depending on the inclination, the disk flux can be a 
factor of two higher or lower compared to a blackbody 
radiating disk. 

— Strong iron line blanketing causes broad (> 100 A) 
spectral features that could be detectable even from 
almost edge-on disks. Disks that are cooler than our 
model (because of a lower mass-accretion rate) could 
even exhibit a spectral break in the optical band due 
to massive line blanketing. 

We stress that these results hold strictly only for our 
particular disk model. In order to arrive at more general 
results a systematic parameter study (disk extent, accre- 
tion rate) of the disk emission is necessary. Also, it needs 
to be investigated in detail how fine the subdivision of 
the disk in a number of rings is necessary in order to 
achieve a computed spectrum with a certain accuracy. 
In addition, deviations from the a-disk model must be 
studied. Another important point will be the inclusion of 
disk irradiation by the X-ray emission from the neutron 
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